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Abstract 

A  phenomenological  characterization  of  hysteresis  in  magnctostrictivc  materials  is  pre¬ 
sented.  Such  hysteresis  is  due  to  both  the  driving  magnetic  fields  and  stress  relations  within 
the  material  and  is  significant  throughout  most  of  the  drive  range  of  magnctostrictivc  trans¬ 
ducers.  An  accurate  characterization  of  the  hysteresis  and  material  nonlincaritics  is  necessary 
to  fully  utilize  the  actuator/sensor  capabilities  of  the  magnctostrictivc  materials.  Such  a 
characterization  is  made  here  in  the  context  of  generalized  Prcisach  operators.  This  yields 
a  framework  amenable  to  proving  the  wcll-poscdncss  of  structural  models  that  incorporate 
the  magnctostrictivc  transducers.  It  also  provides  a  natural  setting  in  which  to  develop  prac¬ 
tical  approximation  techniques.  An  example  illustrating  this  framework  in  the  context  of  a 
Timoshenko  beam  model  is  presented. 
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1  Introduction 

Important  members  in  the  class  of  smart  materials  currently  employed  in  control  applica¬ 
tions  arc  magnctostrictivc  materials.  The  phenomenon  of  magnetostriction  is  defined  as  the 
strain  which  results  when  a  magnetic  material  is  subjected  to  a  magnetic  field.  While  this 
phenomenon  occurs  in  most  ferromagnetic  materials,  only  in  recently  developed  rare-earth 
materials  arc  the  strains  and  forces  sufficiently  large  to  facilitate  their  use  in  actuators  and 
sensors.  This  has  led  to  the  use  of  magnctostrictivc  materials  as  ultrasonic  transducers,  posi¬ 
tioners,  sonar  projectors  (500-2000  IIz)  and  isolators  (5-60  IIz).  They  arc  also  being  considered 
as  actuators  for  controlling  vibrations  in  thick  structures  and  in  heavy  rotating  components 
such  as  milling  machine  bits  (further  details  regarding  applications  in  which  magnctostrictivc 
transducers  arc  employed  as  sensors  and  actuators  can  be  found  in  [8,  16]). 

To  fully  utilize  their  capabilities  as  either  sensors  or  actuators,  the  input  and  output 
characteristics  of  magnctostrictivc  materials  must  be  quantified  in  a  manner  amenable  to 
parameter  estimation  and  control  applications.  As  detailed  in  the  literature  (c.g.,  see  [17,  19]) 
and  outlined  in  the  next  section,  the  input/output  characteristics  of  magnctostrictivc  materials 
arc  inherently  nonlinear  and  display  significant  hysteresis.  To  utilize  the  transducer  responses 
in  a  range  useful  for  full  control  applications,  accurate  mathematical  characterization  of  the 
nonlincaritics  and  hysteresis  is  required. 

The  general  techniques  we  use  here  to  characterize  the  magnetic  nonlincaritics  and  hys¬ 
teresis  arc  based  on  Prcisach-typc  operators.  The  use  of  Prcisach  techniques  for  characterizing 
general  magnetic  hysteresis  is  well-established  in  the  literature  (see  [1,  6,  13,  15,  23,  24,  26] 
and  references  therein)  and  some  aspects  have  been  extended  to  magnctostrictivcs  [2,  9].  In 
[2],  the  classical  Prcisach  model  was  modified  for  Tcrfcnol-D  materials  by  considering  a  char¬ 
acterization  for  the  output  strain  in  terms  of  two  inputs  and  two  Prcisach  kernels.  This  was 
motivated  by  the  observation  that  magnctostrictivc  materials  exhibit  hysteresis  with  respect 
to  both  magnetic  field  and  stress.  Further  coupling  between  the  kernels  and  input  was  included 
in  the  model  proposed  in  [9]. 

In  this  work,  we  consider  the  modeling  of  nonlincaritics  and  hysteresis  in  magnctostrictivc 
transducers  in  the  context  of  the  generalized  Prcisach  or  Krasnosclskii-Pokrovskii  kernels 
developed  in  [6].  Such  generalizations  of  the  classical  Prcisach  theory  arc  motivated  by  the 
goal  of  attaining  a  kernel  which  is  continuous  with  respect  to  both  time  and  shape  parameters 
(as  proven  in  [6],  the  classical  Prcisach  kernel  is  discontinuous  in  both  aspects).  This  provides 
a  framework  in  which  to  prove  wcll-poscdncss  of  models  which  incorporate  magnctostrictivc 
transducers.  Furthermore,  it  yields  a  framework  which  is  natural  for  the  development  of 
approximation  techniques  for  both  simulations  and  parameter  estimation. 

It  should  be  noted  that  the  modeling  of  nonlincaritics  and  hysteresis  through  Prcisach 
techniques  is  phenomenological  rather  than  physics-based.  It  provides  a  mathematical  charac¬ 
terization  of  the  input/ output  relationships  for  magnctostrictivc  materials  in  lieu  of  a  complete 
theory  for  the  electromagnetic  and  magnctomcchanical  properties  of  the  materials.  As  the 
physical  theory  is  advanced,  physics-based  models  may  provide  additional  insights  regarding 
the  actuator/sensor  mechanisms  and  hence  replace  the  Prcisach  models.  Fbr  example,  theory 
based  upon  magnetic  domain  wall  motion,  which  is  used  in  [21,  22]  to  model  hysteresis  in  fer¬ 
romagnetic  materials,  may  be  applicable  to  magnctostrictivc  materials.  At  the  current  time, 
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however,  several  physical  mechanisms  in  magnetostrictive  materials  are  still  not  well  under¬ 
stood,  thus  motivating  the  phenomenological  or  empirical  characterization  of  the  material. 

The  construction  of  magnetostrictive  transducers  and  physical  properties  of  the  magne- 
tostrictivc  materials  are  described  in  Section  2.  The  source  and  form  of  the  nonlinearities 
and  hysteresis  are  detailed  to  illustrate  issues  to  be  addressed  in  the  models.  The  generalized 
Preisach  kernels  and  operators  developed  in  [6]  arc  summarized  in  Section  3.  In  Section  4, 
a  thick  cantilever  beam  with  magnetostrictive  actuators  is  used  as  a  prototype  for  structural 
models  which  incorporate  magnetostrictive  inputs.  Modified  Timoshenko  equations  arc  used 
to  model  the  beam  dynamics  while  generalized  Preisach  or  Krasnosclskii-Pokrovskii  opera¬ 
tors  arc  used  to  characterize  the  inputs.  The  wcll-poscdncss  of  the  model  and  parameter 
estimation  problem  is  established  using  the  framework  of  [4,  6].  Although  the  framework  is 
described  here  in  the  context  of  a  specific,  simple  structure,  it  is  sufficiently  general  to  include 
structures  comprised  of  plates  and  shells  as  well  as  structures  in  which  nonlinear  dynamics 
arc  significant. 

2  Magnetostrictive  Transducers 

To  illustrate  issues  which  must  be  addressed  when  modeling  magnetostrictive  materials,  the 
transducer  depicted  in  Figure  1  is  considered.  This  transducer  is  typical  of  those  currently 
employed  in  applications  and  contains  the  basic  components  required  for  experimental  use  as 
a  control  actuator.  Details  regarding  the  construction  and  performance  of  this  transducer  can 
be  found  in  [17]. 

The  primary  components  of  the  transducer  consist  of  a  magnetostrictive  rod,  a  wound  w'ire 
solenoid,  and  a  cylindrical  permanent  magnet.  In  current  transducers,  the  magnetostrictive 
material  is  typically  composed  of  terbium  and  dysprosium  alloyed  with  iron.  A  commonly 
employed  material  is  Terfenol-D  (Ter:  terbium,  fe:  iron,  nol:  Naval  Ordinance  Laboratory,  D: 
dysprosium)  which  is  constructed  as  a  cylindrical  rod  and  placed  in  the  center  of  the  trans¬ 
ducer.  The  sensor/actuator  capabilities  of  the  material  are  due  to  the  rotation  of  magnetic 
moments  within  the  rod  in  the  presence  of  an  applied  magnetic  field.  In  ferromagnetic  ma¬ 
terials  such  as  Terfenol-D,  moments  are  ordered  and  exhibit  a  high  degree  of  alignment  at 
temperatures  below  the  Curie  point  (regions  in  which  moments  are  highly  aligned  are  termed 
domains).  As  depicted  in  Figure  2a,  the  moments  within  domains  are  primarily  oriented  per¬ 
pendicular  to  the  longitudinal  rod  in  the  absence  of  an  applied  field.  Prestressing  the  rod  with 
the  spring  washer  serves  to  increase  the  number  of  moments  perpendicular  to  the  axis  (see 
Figure  2b)  and  places  the  material  in  compression.  This  latter  objective  is  necessary  due  to 
the  inherent  brittleness  of  Terfcnol-D.  In  the  presence  of  a  magnetic  filed,  the  moments  rotate 
so  as  to  align  with  the  field.  Consequently,  if  the  field  is  applied  in  the  direction  of  the  rod 
axis,  the  moments  rotate  in  the  sense  depicted  in  Figure  2c  and  significant  strains  arc  gener¬ 
ated.  This  is  termed  the  Joule  effect  and  provides  the  actuator  capabilities  of  the  transducer. 
Sensing  is  accomplished  through  the  measurement  of  the  magnetic  fields  which  result  when 
mechanical  stresses  cause  rotations  of  the  moments  within  domains  (Villari  effect).  Details 
regarding  these  effects  can  be  found  in  [19,  25]. 

The  strains  generated  through  an  applied  field  arc  always  positive  since  rotation  of  the 
moments  from  the  prestressed  perpendicular  state  leads  to  an  increase  in  the  rod  length.  As 
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Figure  1.  Cross  section  of  a  typical  Tcrfcnol-D  magnctostrictivc  transducer. 

indicated  in  Figure  3,  the  relationship  between  the  applied  magnetic  flux  or  induction  B  and 
strain  c  is  also  highly  nonlinear  w'ith  saturation  occurring  at  large  field  strengths.  Moreover, 
slight  hysteresis  also  exists  between  B  and  c  at  high  drive  levels  (this  is  not  depicted  in  the 
figure). 

The  generation  of  bidirectional  strains  is  accomplished  through  either  a  DC  current  Iq 
applied  to  the  solenoid  w'hich  surrounds  the  rod,  or  an  enclosing  cylindrical  permanent  magnet 
which  provides  a  biasing  magnetic  induction  Bq.  A  time  varying  current  I{t)  is  then  used  to 
vary  the  induction  in  the  rod  between  0  and  This  provides  the  capability  of  generating 
both  positive  and  negative  strains. 


I  I 
AX2 


(C) 

Figure  2.  Magnetic  domains  in  the  Tcrfcnol-D  rod;  (a)  Orientation  of  moments  within 
domains  in  unstressed  rod  with  no  applied  field;  (b)  Orientation  of  moments  in  prestressed 
rod  with  no  applied  field;  (c)  Orientation  of  moments  in  prestressed  rod  when  field  is  applied 
in  direction  of  longitudinal  rod  axis. 


3 


Figure  3.  Strain  distribution  c  generated  by  an  applied  magnetic  induction  D. 

To  model  the  transducer  for  actuator  and  sensor  purposes,  it  is  necessary  to  characterize 
the  relationship  between  the  current  1  applied  to  the  solenoid,  the  resulting  magnetic  field  //, 
the  associated  magnetic  induction  D  and  finally,  the  generated  strains  c.  Furthermore,  the 
quantification  must  incorporate  the  contributions  due  to  the  permanent  magnet. 

As  detailed  in  [19],  the  magnetic  induction  and  field  are  related  by  the  permeability  fi 
which  is  defined  as 

D 


The  magnetic  field  is  due  both  to  the  solenoid  and  the  permanent  magnet.  The  magnetic 
induction  from  the  permanent  magnet  is  approximated  by  Dq  =  /.iIIq  while  Ampere’s  law 
yields  B  =  finl^  where  n  is  the  number  of  turns  per  unit  length  in  the  solenoid,  as  the 
approximate  magnetic  induction  due  to  the  solenoid  (these  arc  approximate  since  edge  effects, 
air  gaps,  etc.,  arc  neglected). 

For  Tcrfcnol-D,  the  permeability  /.i  is  highly  nonlinear  and  exhibits  significant  hysteresis 
as  indicated  by  the  induct ion/ficld  relations  depicted  in  Figure  4a.  As  discussed  in  [2],  this 
hysteresis  is  manifested  wdth  respect  to  both  the  applied  magnetic  field  and  stresses  within 
the  magnctostrictivc  material.  To  indicate  the  latter  contribution,  the  permeability  in  mag- 
nctostrictivc  applications  is  often  denoted  by  ,  The  hystcrctic  relationship  between  the 
magnetic  field,  magnetic  induction  and  material  stress  arc  then  inherently  manifested  in  the 
field-strain  relations  as  shown  in  Figure  4b. 

As  indicated  by  the  preceding  discussion,  the  nonlinear  relationships  bctw'ccn  the  applied 
current  and  magnetic  induction  arc  augmented  by  nonlincaritics  in  the  behavior  of  the  magnc¬ 
tostrictivc  materials.  Fbr  example,  experimental  results  in  [11,  14]  indicate  that  the  Young’s 
modulus  for  Tcrfcnol  is  dependent  upon  the  applied  magnetic  field  w'hich  partially  accounts 
for  the  dependence  of  magnetic  hysteresis  on  the  material  stress.  Furthermore,  experimen¬ 
tal  results  in  [11,  14]  demonstrate  that  other  material  properties  such  as  magnet omcchanical 
coupling  coefficients  arc  highly  sensitive  to  operating  conditions  such  as  prestress  level,  AC 
drive  levels,  operating  frequencies  and  temperature.  To  provide  a  framework  amenable  to 
characterization  of  these  relationships,  we  summarize  next  necessary  theory  regarding  gener¬ 
alized  Prcisach  operators.  The  application  of  this  theory  to  magnctostrictivc  materials  in  a 
structural  setting  is  considered  in  the  final  section. 
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(a)  (b) 


Figure  4.  (a)  Relationship  between  the  magnetic  field  strength  II  and  the  magnetic  flux 
density  B:  (b)  Applied  magnetic  field  II  and  resulting  strain  distribution  c. 


3  Preisach  and  Krasnoselskii-Pokrovskii  Operators 


In  this  section,  we  provide  basic  theory  regarding  Preisach- type  operators  in  a  Hilbert  space 
setting.  Classical  Preisach  kernels  and  operators  are  defined  first  and  are  then  extended  to 
kernels  of  Krasnoselskii-Pokrovskii  type  to  attain  desired  continuity  properties.  The  material 
in  this  section  summarizes  theory  from  [6,  7],  and  the  reader  is  referred  to  those  references  for 
further  details. 

To  motivate  the  general  kernels  used  later,  we  first  illustrate  with  a  single  delayed,  relay 
operator  k.  This  kernel  is  characterized  in  terms  of  crossing  times  T{t)  defined  by 


r(<)  =  {t)€  (0,  T]  I  u(>7)  =  Sj  or  u(rj)  =  sj) 


where  s  =  (51,52)  are  points  in  the  Preisach  half  plane 

5  =  ^5  €  IR^  I  5  -  (51, 52),  5i  <  52I 

and  u  denotes  an  input  function.  The  values  51,52  are  threshold  values  for  the  multivalued 
kernel  as  reflected  in  the  definition 


[k{u,om 


[^’5(1^01(0)  if 

-1  if 

+1  if 


r(0  =  0 

r(i)  ^  0  and  w(maxr(i))  =  5i 
T(i)  ^  0  and  w(maxr(^))  =  52  . 


A  depiction  of  this  kernel  is  given  in  Figure  5.  The  starting  value 


[i(u,e)](o) 


“1  if  14(0)  <  5i 

^  if  5i  <  u(0)  <  52 
+  1  if  14(0)  >  52 


defines  the  initial  state  of  the  kernel  in  terms  of  the  parameter  {  €  {—1,1}. 
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The  output  remains  on  a  branch  until  a  threshold  is  reached  in  the  monotonically  increasing 
input  u.  At  that  point,  the  output  jumps  to  the  other  saturation  value  and  remains  there 
until  the  other  threshold  value  is  reached.  For  example,  an  output  response  starting  with  a 
value  of  —1  will  retain  that  value  until  u{t)  reaches  52.  The  output  then  jumps  to  -fl  until 
the  threshold  value  of  5 1  is  reached. 

The  classical  Preisach  operators  are  then  defined  in  terms  of  parallel  collections  of  these 
single  relay  operators.  To  this  end,  we  let  M  denote  the  set  of  all  finite,  signed  Borel  measures 
on  S  and  let  /  be  a  Borel  measurable  function  mapping  S  ->  {—1,1}.  For  u  €  C[d^T]  and 
€  My  the  Preisach  operator  is  defined  by 

[>.(«,/)](<)  = 

The  goal  in  the  parameter  identification  problem  is  to  estimate  the  measure  1/  so  that  the 
model  response  “fits”  experimental  data  in  a  least  squares  sense. 

While  this  provides  an  operator  which  is  useful  for  many  applications,  this  classical  defi¬ 
nition  does  not  yield  a  kernel,  and  hence  operator,  which  is  continuous  with  respect  to  either 
time  or  parameters.  Specifically,  as  proven  in  [6],  the  mapping  in  time 

[kiu,()]{t) 


and  the  parameter  space  mapping 

s  1-^  [ks{u,  0](<) 

are  discontinuous  for  the  classical  Preisach  kernel  ks(uyC)'  Continuity  in  time  is  important 
from  a  physical  perspective  while  continuous  parameter  dependence  is  crucial  for  the  devel¬ 
opment  of  practical  parameter  estimation  methods. 


Figure  5.  Single  Preisach  relay  operator  wdth  threshold  values  51,52. 
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To  avoid  the  difBcultics  associated  with  the  discontinuous  mappings,  a  Krasnosclskii- 
Pokrovskii  kernel  of  the  type  discussed  in  [6]  is  employed.  This  kernel  is  somewhat  less 
general  than  the  influence  operators  considered  in  [23]  and  arises  as  an  extension  of  smoothed 
Preisach  operators.  These  operators  differ  from  the  previously-defined  Preisach  operator  in 
the  manner  through  which  an  envelope  of  admissible  paths  is  defined.  In  this  case,  an  envelope 
is  provided  by  translates 

r,!  =  r(x  -  s,) 

=  r{x  -  S2) 

of  a  Lipschitz  continuous  ridge  function  r(a;)  as  depicted  in  Figure  6.  For  monotone  inputs 
Ujn  €  CfO,!"],  a  monotone  output  operator  is  defined  by 


max{^,  r(u,n(0  “  •52)}  if  Um  is  non- decreasing 
min{^,r(u„i(f)  —  51))  if  is  non-increasing . 


In  terms  of  this  operator,  a  kernel  is  defined  for  piecewise  monotone  inputs  Upy^  €  C'[0,T]  fl 
5i,j[0,T],  where  5i,j[0,T]  is  the  set  of  piecewise  linear  splines  with  j  knots  in  [0,T],  in  the 
following  manner.  The  initial  value  of  the  operator  is  taken  to  be  7?o  =  {.  A  kernel  ks  is  then 
defined  recursively  on  each  subinterval  by 


[ks(^Upjfly  o]it)  =  [n  )](<).  [4-1, f.l  (3.1) 

where  71^  =  ^  =  f?*  *  *  The  input  and  action  of  this  kernel  are  illustrated 

in  Figure  7.  This  provides  a  definition  of  the  kernel  useful  for  computational  algorithms  in 
which  inputs  arc  discretized  in  terms  of  a  piecewise  linear  basis.  This  definition  is  readily 
extended  to  arbitrary  u  €  C[0,T]  through  standard  density  arguments  as  detailed  in  [23]. 

It  is  natural  in  the  applications  considered  here  to  formulate  the  hystcrctic  input  operator 
in  a  manner  commensurate  with  a  weak  model  formulation.  We  thus  consider  a  space  of  test 
functions  V"  and  state  space  II  in  w'hich  V'  is  continuously  and  densely  embedded  in  II  and 
forms  a  Gclfand  triple 

V  -  ir  ^ 


Figure  6,  Hysteresis  envelope  provided  by  the  translates  ri  and  r2  of  the  ridge  function  r(rr). 
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Figure  7.  (a)  Picccwisc  monotone  input;  (b)  Output  from  the  Krasnoselskii-Pokrovskii  kernel 
in  response  to  a  piecewise  monotone  input. 


The  following  theorem  from  [6]  summarizes  the  continuity  properties  of  the  kernel  and 
quantifies  the  resulting  input  operator  Bu- 

Theorem  1.  Let  kg  denote  the  kernel  defined  in  (3.1)  and  let  S  denote  the  closure  of  S  in  IR^. 

(1)  For  each  s  £  S  and  (  £  {—1, 1),  the  kernel  satisfies  ks{'y  (^)  :  f7[0,  T]  C'[D,  T] . 

(2)  For  each  u  £  6'[0,T],^  £  (—1, 1}  and  t  £  C[0,T],  the  map  s  h->  [A:5(u,  (^)](i)  is  continuous 

from  S  to  IR. 

(3)  For  g  £  V*  and  Jy'  £  the  control  influence  operator  i?,,  defined  by 

[Z?.(«,0](‘)  =  =  9 1  [hM]it)d>y{s)  (3.2) 

satisfies  £  L^((0,  T);  V"*). 


Remark  1.  The  kernel  kg  and  control  influence  operator  R*,  defined  in  this  manner  can  be 
readily  incorporated  in  a  framework  for  proving  model  w'ell-posedness  and  the  existence  of 
a  measure  which  minimizes  an  appropriate  parameter  estimation  functional.  Moreover,  as 
illustrated  in  Figure  7,  this  kernel  yields  a  technique  for  characterizing  the  nested  hysteresis 
curves  which  arc  common  in  magnetostrictive  applications. 


4  Characterization  of  Hysteresis  for  Structural  Models 
with  Magnetostrictive  Transducers 

To  illustrate  the  use  of  the  operator  framew'ork  described  in  the  last  section  to  characterize 
hysteresis  in  magnetostrictive  materials,  we  consider  the  modeling  of  magnetostrictive  trans¬ 
ducers  mounted  to  a  thick  cantilever  beam  as  depicted  in  Figure  8.  The  transducers  arc 
considered  to  be  mounted  to  the  clamps  at  the  fixed  edge  of  the  beam  so  that  mass  load¬ 
ing  from  the  actuators  themselves  is  minimized.  A  rigid  bar  is  used  connect  the  end  of  the 
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Figure  8.  Cantilever  beam  with  magnetostrictive  actuators. 


Terfenol-D  rod  in  the  transducer  to  the  beam.  By  driving  the  transducers  out-of-phase,  bend¬ 
ing  moments  are  generated  in  a  manner  which  can  be  used  to  attenuate  beam  vibrations.  As 
described  in  [12],  this  experimental  setup  has  been  used  in  initial  experiments  to  determine 
the  potential  of  magnetostrictive  transducers  as  structural  actuators.  Due  to  limitations  in 
models  and  control  laws,  driving  currents  in  the  experiments  were  restricted  to  a  range  in 
which  linearized  results  could  be  employed.  Even  in  this  restricted  regime,  the  results  of  [12] 
demonstrate  the  utility  of  the  magnetostrictive  transducers  for  structural  applications. 

It  should  be  noted  that  a  beam  is  considered  here  both  due  to  its  previous  experimental  use 
and  the  relative  simplicity  of  the  resulting  model.  The  operator  techniques  are  quite  general, 
however,  and  can  be  applied  in  a  similar  manner  to  structures  comprised  of  plates  or  shells  as 
well  as  structures  undergoing  large  deformations  which  leads  to  nonlinear  models  (e.g.,  von 
Karman  models). 

For  modeling  purposes,  we  take  the  beam  to  have  length  width  b  and  thickness  h.  The 
density,  Young’s  modulus,  Kelvin- Voigt  damping  coefficient  and  air  damping  coefficient  for  the 
beam  are  denoted  by  and  7,  respectively.  The  cross-sectional  area  of  the  Terfenol 

rod  is  denoted  by  wffiilc  the  Young’s  modulus  and  damping  coefficient  for  the  Terfenol 

rod  are  denoted  by  cf,.  The  length  and  width  of  the  connecting  bar  arc  denoted  by  ^ 

and  br,  respectively,  while  the  bar  density  is  given  by  p,.. 


4.1  Strong  Form  of  Beam  Model 

The  Timoshenko  equations 


(t,x} 


0<x<£, t>0 


(4.1) 


w(t,0)  =  a(t^O)  =  0  1 

A^(t,£)  =  Q(t,£)  =  0  J  ’ 


where  w  and  a  denote  the  transverse  displacement  and  cross-sectional  rotation,  respectively, 
arc  used  to  model  the  beam  dynamics.  Exogenous  surface  forces  to  the  beam  arc  denoted  by 


f{tyx)  while  p{x)  is  the  composite  density  of  the  structure.  The  cross-sectional  area  of  the 
beam  is  given  by  A(x)^  while  r{x)  =  I{x)/A{x)  where  I(x)  is  the  moment  of  inertia  of  the  cross- 
sectional  area.  Note  that  these  quantities  are  spatially  variable  due  to  the  nonhomogcncity 
in  the  region  of  the  connection  rod.  Finally,  the  internal  bending  moment  and  shear  force  are 
given  by  M{t^x)  and  Q{tyx)^  respectively,  while  the  external  bending  moment  generated  by 
the  magnctostrictive  actuators  is  denoted  by  Mmagittx). 

To  determine  appropriate  functional  forms  for  the  density,  internal  moment  and  shear  force, 
the  structural  contributions  due  to  the  connecting  bar  and  Terfenol  rod  must  be  quantified. 
VVe  assume  here  that  the  connecting  bars  arc  perfectly  rigid  and  contribute  mass  to  the  beam 
but  do  not  affect  the  bending  moments  (we  neglect  air  resistance  to  the  bars).  The  actuator 
and  Terfenol  rod  arc  considered  to  be  supported  from  the  boundary  clamps  so  they  do  not 
contribute  mass  to  the  beam.  The  Terfenol  rod  is  assumed  to  contribute  an  clastic  stress 
w'hich  is  uniform  across  the  cross-sectional  area  of  the  rod. 

Under  the  assumption  of  uniform  cross-sectional  strains  in  the  magnctostrictive  rods,  the 
density,  stiffness  and  Kelvin- Voigt  damping  parameters  for  the  structure  arc  then  taken  to  be 


p(x)  =  pihb  +  2p,.h4rXrod{x) 

EI{x)  =  ^  {h/‘2  +  i,.f  xM  (4.2) 

cul{x)  =  +  24„.,c"  ih/2  +  £,f  xM 

where  the  location  of  the  rods  is  delineated  by  the  characteristic  function  Xrod  which  has  a 
value  of  1  in  the  region  covered  by  the  connection  bar  and  is  0  elsewhere.  The  internal  moment 
and  shear  arc  then  given  by 


M{t,x)  =  El{x)^{t,x)  +  culix)-^^{t,x) 


Q{ty  a:)  =  KAG{x)3{ty  x)  -f  «/lcg(a:)-^(t,  a:) 


(4,3) 


where  /c  is  a  correction  factor  which  accounts  for  the  fact  that  the  outer  surface  of  the  beam 
cannot  support  a  shear  stress,  ^(a:)  is  the  shear  modulus  and  cg(a:)  represents  resistance  to 
the  shear  strain  rate.  Finally,  the  shear  deformations  arc  defined  by 


0{tyX) 


div 

dx 


{tyX)-\-a(tyX), 


It  should  be  noted  that  the  contributions  due  to  the  connection  rods  arc  dependent  upon 
the  exact  experimental  setup  and  different  assumptions  and  models  can  also  be  used  to  in¬ 
corporate  the  passive  rod  contributions.  In  all  cases,  the  picccwdsc  constant  parameters 
p{x)yr{x)y  EI{x)yG{x)yCuI{x)yCQ{x)  cLiid  constant  parameter  7  must  be  estimated  through 
a  least  squares  fit  to  experimental  data  to  attain  a  reasonable  model  for  the  specific  experi¬ 
mental  device  (The  values  determined  by  (4.2)  cannot  be  used  with  certainty  when  modeling 
the  experimental  apparatus  due  to  inaccuracies  in  manufacturer  specifications,  etc.;  however. 
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they  can  be  used  as  initial  values  for  the  optimization  routine.)  Note  that  in  the  moment  and 
shear  expressions  (4.3)  and  equations  (4.1),  parameters  have  been  combined  to  yield  a  single 
value  to  be  estimated. 

To  characterize  the  external  moment  generated  by  the  magnetostrictive  transducers,  we 
w'ill  assume  that  an  offset  DC  current  Iq  is  used  to  provide  the  magnetic  field  bias  necessary 
to  attain  bidirectional  strains.  Ampere’s  law  then  yields  the  magnetic  induction 

B{t)  ~  +  lo]  (4.4) 

where  again,  n  is  the  number  of  terms  per  unit  length  in  the  solenoid,  I{t)  is  the  current 
applied  to  the  solenoid,  and  is  the  permeability.  As  noted  in  Section  2,  the  permeability 
exhibits  both  stress-dependent  nonlinearities  and  hysteresis  due  to  the  magnetic  field. 

The  strain  can  be  related  to  the  magnetic  induction  through  the  nonlinear  equation 

c{t)  =  fCBBit)  (4.5) 

w'hcre  JCb  depends  upon  the  magnetic  induction  D,  Moreover,  the  external  stress- strain 
relation 

a{t)  -  E^c(t)  (4.6) 

is  also  nonlinear  since  the  Young’s  modulus  for  the  Terfcnol-D  rod  can  be  highly  dependent 
upon  the  applied  magnetic  field.  Finally,  the  external  moment  generated  by  the  Terfenol-D 
rod  is 

M,nug{t,  X)  =  2A,nag{ir  +  hl2)(j{t)Xrod{x)  (4.7) 

where  tr  denotes  the  length  of  the  connection  bar  and  A, nag  is  the  cross-sectional  area  of  the 
Terfenol  rod. 

Combination  of  (4.7)  with  (4.4)-(4.6)  will  yield  a  relationship  for  the  external  moment  in 
terms  of  the  applied  current  l{i).  However,  this  relationship  is  not  directly  useful  for  modeling 
transducer  dynamics  since  the  parameters  and  E^  are  highly  nonlinear  and  exhibit 

significant  hysteresis. 

An  alternative  is  to  characterize  the  external  moment  through  the  fitting  of  hysteresis  and 

material  nonlincaritics  in  terms  of  the  Preisach  kernels  described  in  Section  3.  Specifically, 

the  external  moment  generated  by  the  magnetostrictive  transducer  can  be  described  by  the 
relation 

M,nagit,  a)  =  2Anag{^r  d-  h/2)[P^{I,  0]{t)Xrodix)  (4.8) 

where 

[PM,0]it)  =  l[ksiI,0mMs)-  (4.9) 

The  kernel  ks  is  the  extension  of  that  defined  in  (3.1)  to  I  £  C[0,!r].  In  essence,  the  method 
characterizes  the  hysteresis  and  material  nonlinearitics  through  curve  fitting  in  terms  of  a 
measure  u  which,  for  a  given  set  of  experimental  operating  conditions,  is  estimated  through 
a  least  squares  fit  to  data.  Due  to  its  generality,  the  technique  can  be  used  to  characterize 
a  wide  variety  of  transducer  responses.  Furthermore,  the  technique  provides  a  framework 
suitable  for  parameter  estimation  and  subsequent  controller  design. 
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4.2  Weak  Form  of  Beam  Model 

As  noted  when  (4.3)  and  (4.8)  arc  employed  in  (4.1),  the  use  of  the  strong  form  of  the  beam 
model  leads  to  the  differentiation  of  discontinuous  material  parameters  and  inputs.  It  also 
necessitates  the  use  of  high-order  approximating  elements.  To  alleviate  these  difficulties  and 
provide  a  framework  amenable  to  analysis  and  approximation,  we  consider  a  corresponding 
weak  form  of  the  model. 

For  this  system,  the  state  is  taken  to  he  y  =  {w,  q)  in  the  Hilbert  space  II  =  L^(0,^)  x 
L^(0,  t)  with  the  inner  product 

(4>,\P')^=  f  pSitl^idx-h  f  pr(f)2i^2dx 

Jo  Jo 

where  4>  —  (<^1,(^2), '5  =  (V'i)V'2)-  The  space  of  test  functions  is  V"  =  //|(0,  f)  X  //^(O,  £), 
7/2(0,^)  =  [(f)  ^  77^  (0,7)  I  <^(0)  =  0),  with  the  inner  product 

(4>,  =  f\6[  -1-  ^2)(V^;  +  xl^2)  dx  +  dx  . 

Jo  Jo 

It  should  be  noted  that  with  these  choices,  V"  is  continuously  and  densely  embedded  in  77  and 
that  and  77  form  a  Gelfand  triple;  that  is 


V  77  -  ir  -4  V* . 


A  weak  form  of  the  modeling  equations  is  then 

pw^)\  dx  -j-  ^;wxb\  dx  +  J  Qi!^)[  dx  ^  J  prail)2dx  -|-  J  +  J  Q^-2  dx 


f 


r 

/  ^Imagy-‘2^^  T  /  fVj  dx 
Jo 


(4.10) 


for  all  'P  =  (V’i?V'2)  €  V^  In  this  form,  derivatives  are  transferred  onto  suitably  smooth 
test  functions.  This  alleviates  the  difficulties  associated  with  the  discontinuities  and  reduces 
smoothness  requirements  on  approximate  solutions. 


4.3  Model  Well-Posedness 

To  provide  a  framework  amenable  to  proving  the  well-posedness  of  the  model  and  parameter 
estimation  problem,  it  is  advantageous  to  formulate  the  model  in  terms  of  sesquilinear  forms 
and  the  operators  which  they  generate.  To  this  end,  recall  that  y{t)  ~  (u;(<,  •),  a(i,  •))  and 
define 


c^i(^y)(y(0)'‘^)  =  f  EIa'il>2dx+  [  kAG{io' -h  -h  11^2)  dx 
Jo  Jo 

^2iQ){y{t),  'P)  =  f  coloiiy^  dx-h  f  kAcq{w'  A  cv)(t/-’(  +  ^-2)  dx  -j-  f  'ywrfyi  dx 


(4.11) 
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where  q  —  (p,  r,  EI^  G,  c^)/,  cq^j)  is  considered  in  an  admissible  parameter  space  Q,  It  can  be 
directly  verified  that  ai  and  cr2  satisfy  the  boundedness  and  ellipticity  conditions 


(III) 

ki(9)(^,'P)l<c,|$|i/|'P|v'  , 

Cl  ^  IR 

(112) 

Rcai(9)($,§)>C2|$|2  ,  C2 

>  0 

(113) 

M9)(^,'P)|<C3l€)K'hPk  , 

C3  e  E 

(114) 

Ilccr2(9)($,§)  +  C4|$|^  >  Co 

ml  ,  C4 

for  all  €  V"  where  the  constants  ci,  •  •  • ,  C5  arc  independent  of  the  parameters  q  €  Q. 

From  (4.8),  it  follows  that  the  transducer  contributions  can  be  represented  in  terms  of  the 
operator  :  U  — >  V'*  by 


Om,  =  ([-P.(m,  e)](%, 

-  [^.(w.  01(0  r  +  hl2)i^Xroi  dx 

Jo 

where  U  =  ITl,  (*,*)v/*  is  the  usual  duality  product  and  P,/  is  defined  in  (4.9).  Note  that 
fromJThcorcm  1,  it  follows  that  Bu{u,^)  G  L^((0,  T);  V^*).  Finally,  with  the  definition  f{t)  = 
(p“V(i,  •),0),  we  can  write  the  weak  form  (4.10)  in  the  abstract  variational  form 

{Ht)y  '5)1'*^  +  <^2{q)iy(t),  'P)  +  «^i(9)(y(<),  'P)  =  {[B„{u,0W)  +  fi*),  '^)v\v 

y(0)  =  2/0  ,  i'(O)  =  yi 


for  all  'P  G  V. 

An  equivalent  system  can  be  obtained  by  invoking  the  boundedness  of  cri,cr2  to  define 
operators  /1,(^)  G  C{\'\V*)  =  1,2,  by 


vp)j,,.  j,,  =  <T2(9)($,  vP)  ’  ^  ^  • 

In  operator  form,  the  equations  governing  the  beam  dynamics  then  have  the  form 


y{t)  +  A2{q)y{t)  +  Ai{q)y{t)  =  lD^(u,  0](<)  +  /(<)  (4.13) 

in  y*. 

The  subsequent  result  concerning  the  existence,  uniqueness  and  regularity  of  solutions 
follows  directly  from  Theorem  2.1  and  Remark  2.1  of  [4]. 


Theorem  2.  Let  Q  denote  a  compact  metric  space  and  j\4  denote  the  set  of  finite,  signed 
Borcl  measures  on  5.  Consider  inputs  u  G  G[0,r]  and  /  G  {—1,1}).  Finally,  let  cri,cr2 
be  given  by  (4.11)  and  hence  satisfy  (II1)-(II4).  For  each  (5,  j^)  G  Q  x  yV4,  there  then  exists  a 
unique  solution  y  to  (4.12)  which  satisfies 

2/eC((0,r);lOcI^((0,T);V') 

ye(7((0,r);//)nL^((0,T);y) 

j/eL^((0,T);y*). 
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Remark  2.  As  noted  in  Remark  4.2  of  [8],  the  regularity  results  in  Theorem  2  can  actually 
be  strengthened  to  vield 

y€C((0,r);lO 

yecmryj!). 


4.4  Parameter  Estimation 


The  goal  in  the  parameter  estimation  problem  is  to  determine  material  parameters  q  — 
(p,  r,  Zi?/,  (7,  Cx)/,  cg,  7)  and  a  measure  u  given  data  measurements  z  from  some  observable 
subspace  Z  of  the  state  space.  The  form  of  the  quadratic  functional  to  be  minimized  depends 
on  the  experimental  data  which  is  available.  For  time  domain  data  consisting  of  position, 
velocity  or  acceleration  measurements  at  points  x  on  the  beam,  an  appropriate  functional  is 


=  13 


(4.14) 


With  s  =  0, 1  or  2.  Other  functionals  incorporating  cither  time  or  frequency  domain  data  can 
be  found  in  [8].  In  each  case,  the  minimization  is  performed  subject  to  y  satisfying  (4.12)  for 
q  £  Q  and  ly  G  A  C  M  where  A  is  an  appropriate  class  of  measures. 

As  noted  previously,  the  parameters  p,  r,  EI^  (7,  cdI^  cq  arc  assumed  to  be  picccwdsc  con- 
stant  with  partition  points  at  the  connection  rod  edges.  The  air  constant  7  is  constant,  and 
all  seven  parameters  arc  positive.  The  admissible  parameter  space  Q  is  then  taken  to  be  a 
compact  subset  of  the  metric  space  Q  =  x  IR  with  elements  piecewise  constant 

bctw'ccn  partition  points  and  satisfying  the  positivity  constraints. 

An  appropriate  choice  for  v4  is  a  set  of  probability  measures  with  a  metric  which  yields 
convergence  in  distribution.  Specifically,  let  Sa  be  a  compact  subset  of  §  and  let  V{Sa) 
denote  the  set  of  Borcl  probability  measures  u  on  5a.  The  space  A  is  then  defined  to  be 
'P(5a)  endowed  with  the  Prohorov  metric  w'hich  is  defined  for  ui^v'2  €  'P(5a)  by 


=  inf  {e  >  0  1  jyi{F)  <  iy2{F^)  s,  F  closed,  F  C  5a} 

(sec  [10]  for  details  regarding  this  metric).  Here  F^  denotes  an  e  neighborhood  for  F. 


Remark  3.  As  detailed  in  [3,  6,  7,  10],  the  space  A  satisfies  the  following  properties. 

(a)  Let  C(5a)  denote  the  space  of  continuous  functions  on  5a.  Since  P(5a)  C  [(7(5a)]*, 
convergence  in  the  Prohorov  metric  is  equivalent  to  weak*  convergence;  that  is 

74  1/  f  fdi/f,  ->  f  fdi/ 

for  /  €  (7(5a). 

(b)  ^  is  a  compact  metric  space. 
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(c)  Consider  the  control  influence  operator  defined  in  (3.2)  and  let  y{q,u)  be  a  solution  to 
(4.12).  For  the  functional  J{q,v)  defined  in  (4.14)  with  fixed  q^Q,  the  map  u  -4  J{q,u) 
is  weak*  lower  semi  continuous  from  P(5a)  to  111. 

The  following  theorem  taken  from  [6]  specifics  conditions  under  which  the  parameter  esti¬ 
mation  problem  with  the  magnctostrictivc  material  inputs  is  well-posed.  When  combined  with 
theory  from  [5],  it  provides  a  framework  for  numerically  estimating  physical  and  hysteresis 
shape  parameters  through  a  least  squares  fit  to  experimental  data. 

Theorem  3.  Let  denote  the  magnctostrictivc  input  operator  defined  in  (3.2)  and  y{q^u) 
be  the  solution  to  (4.12).  If  5a  is  a  compact  subset  of  5,  then  there  is  a  probability  measure 
€  'P(5a)  which  solves  the  minimization  problem 

J{q,PQ)  =  inf  J{q,u) 

veV{Sj^) 

for  the  functional  J(g, //)  defined  in  (4.14). 


5  Concluding  Remarks 

This  paper  addresses  the  characterization  of  material  nonlincaritics  and  hysteresis  inherent  to 
magnctostrictivc  materials  at  middle  to  high  range  drive  levels.  The  hysteresis  is  induced  by 
both  magnetic  fields  and  stresses  within  the  material  while  additional  nonlincaritics  arise  in 
the  strain-magnetic  induction  relations  and  stress-strain  relations.  Some  of  the  mechanisms 
leading  to  these  nonlinear  input/output  relations  have  been  modeled  in  terms  of  electromag¬ 
netic  and  magnctomcchanical  theories  [20,  21,  22].  However,  several  mechanisms  governing 
magnctostrictivc  properties  arc  still  not  completely  characterized  by  physics-based  models. 
This  motivates  a  mathematical  characterization  in  terms  of  phenomenological  Prcisach  tech¬ 
niques. 

Prcisach  models  arc  empirical  in  the  sense  that  they  can  be  used  to  mathematically  rep¬ 
resent  hysteresis  curves  in  terms  of  shape  parameters  determined  through  least  squares  fits 
to  experimental  data.  Since  the  models  arc  not  directly  derived  from  the  physics  of  the  sys¬ 
tem,  they  can  be  used  to  model  dynamics  in  which  the  underlying  physics  is  not  thoroughly 
understood.  The  price  paid  for  this  generality  is  the  loss  of  insight  which  can  sometimes  be 
provided  by  a  physic-based  model. 

In  this  work,  Krasnosclskii-Pokrovskii  kernels  arc  used  to  characterize  magnctostrictivc 
inputs  to  a  flexible  structure.  As  detailed  in  [6],  such  kernels  arc  advantageous  over  classical 
Prcisach  kernels  due  to  the  property  that  they  arc  continuous  with  respect  to  both  time 
and  shape  parameters.  The  operators  generated  via  these  kernels  arc  then  used  to  formulate 
the  models  in  a  Hilbert  space  framework.  Wcll-poscdncss  of  the  models  is  obtained  through 
the  theory  of  [4]  while  existence  of  optimal  physical  parameters  and  an  optimal  measure  for 
characterizing  the  hysteresis  curve  follows  from  the  theory  of  [3,  5,  6].  Hence  for  structural 
systems  with  magnctostrictivc  transducers,  this  characterization  technique  provides  a  model 
which  is  amenable  to  analysis  and  approximation.  In  future  investigations,  the  practical 
efficacy  of  this  model  w'ill  be  tested  through  validation  experiments. 
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